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THE DISTRIBUTION OF PHASE SHIFTS FOR SEMICLASSICAL 
POTENTIALS WITH POLYNOMIAL DECAY 

JESSE GELL-REDMAN AND ANDREW HASSELL 


Abstract. This is the third paper in a series OE] analyzing the asymptotic distribution 
of the phase shifts in the semiclassical limit. We analyze the distribution of phase shifts, 
or equivalently, eigenvalues of the scattering matrix, Sh(E), for semiclassical Schrodinger 
operators on which are perturbations of the free Hamiltonian by a potential V with 
polynomial decay. Our assumption is that V{x) ~ as a: ^ oo, for some a > d, 

with corresponding derivative estimates. In the semiclassical limit h ^ 0, we show that 
the atomic measure on the unit circle defined by these eigenvalues, after suitable scaling 
in h, tends to a measure /i on Moreover, is the pushforward from R to R/27rZ = 
of a homogeneous distribution u of order /3 depending on the dimension d and the rate 
of decay a of the potential function. As a corollary we obtain an asymptotic formula for 
the accumulation of phase shifts in a sector of 

The proof relies on an extension of results in m on the classical Hamiltonian dy¬ 
namics and semiclassical Poisson operator to the class of potentials under consideration 
here. 


1. Introduction 

Consider a semiclassical Schrodinger operator 

Hh := + V -E 

on where A = — is the positive Laplacian, E is a real constant and V : R 

is a smooth real-valued function satisfying 

( 1 . 1 ) V{x) = + W{x), X G R'^, X = Ai where |lU(x)| = 0 (|x|“^“''''^^), 

|x|“ |x| 

for |x| large, some a > 1 and some e > 0. (For our main theorem, we will require a > d, and 
W to satisfy ‘symbolic’ derivative estimates as in ( 113 , but for some of our intermediate 
results a > 1 and dm will be sufficient.) 

Under these assumptions, the (relative) scattering matrix Sh exists and is a unitary 
operator on L^(S‘^“^); Sh is given on ^ G C'‘’”(S‘^“^) as Shtj) = where '0 is the 

unique function such that there is a solution to HhU^i, = 0 satisfying 

( 1 . 2 ) r = |x|. 

The difference Sh — Id is a compact operator on and thus the spectrum of Sh lies 

on the unit circle, is discrete, and accumulates only at 1. Setting 7 = (d — l)/(a — 1), we 
define the (infinite) atomic measure ^h on the circle which acts on / G \ 1) by 

(1.3) = E /(e^*^-'®) 

^spec{Sh) 

for some enumeration of the eigenvalues of Sh, repeated according to their multiplic¬ 

ity. The Pn,h € [OjTt) ^re called the ‘phase shifts’ of Hh- 
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Main Theorem: Let I be any open interval ofE>^ containing 1. Assume that V satisfies 
(HH) for some a > d, and that W in (HU satisfies the additional derivative estimates 

(1.4) =0(|xr(“+l'=l+'=)) V fee |x|^oo. 

Then the measures pLh converge in the weak-* topology onEf\I to a measure pL. Moreover, 
pL is the pushforward via the map M —> '^jTn'L ~ of a homogeneous measure 

_ I (l/(27r))'^-iai6l-^+i for 9 > 0 
~ { {l/{2TT)y-^a2\9\-'^+^ for9<0. 

Concretely, for every f G C'comp(^^ \ -^): have 

\imUh,f)=f /(e*®)d/J. 

h^O Jgi 

The constants ai, a 2 are given in (14.111) below. 

The Main Theorem is proven at the end of Section 0 modulo the proofs of subsequent 
technical lemmas. 

It follows from the Main Theorem that eigenvalues of Sh accumulate in sectors of the 
unit circle at a rate of Indeed, defining a sector on the circle by choosing 

angles 0 < (/)o < </i’i < 27r, and letting 

N{(j)Q,(l)i) = #{n : 4>q < f3h,n < 4>i mod 27r}, 
the main theorem implies the following. 

Corollary 1.1. Assumptions as in the Main Theorem, the number of eigenvalues in a sector 
satisfies 

Nicfo, f " dpi)il + o(l)). 

J 

This result can be taken as an analogue of the Weyl asymptotic formula, reviewed below 
in Section 15.11 It is also proven at the end of Section O 

The Main Theorem is proven, following [53], via analysis of the traces of the operators 
— Id. The fact that these operators are trace class is shown for example in [22] ; in Section 
m below, we prove a precise asymptotic formula for the trace which gives its leading order 
behavior in h as —>■ 0. The trace of — Id is equal to h~°‘'^{pLh,Pk{z)) with Pk{z) = — 

where pLh is the measure defined in the Main Theorem, and our asymptotic formula for the 
trace of Sj^ — Id shows that the Main Theorem holds for these special values of /. Note 
that pk(z) does not, strictly speaking, satisfy the assumptions of the Main Theorem, as 
its support contains 1; in fact, we prove that the conclusion of the theorem holds on the 
Banach space of continuous functions vanishing to first order at I — see Sections [3] and [Sj 
To conclude that the Main Theorem holds we show in Section [5] that the measures pLh are 
continuous on this space of continuous functions, which contain the span of the pk as a dense 
subset, and use an approximation argument to obtain the formula in the Main Theorem. 

The trace of — Id is obtained via analysis of the Schwartz kernel of Sh and its powers. 
By |12j . with previous results for example in [niiniiis], the operator Sh is a semiclassical 
Fourier Integral Operator whose canonical transformation is the total sojourn relation, 
which is a map from incoming rays to outgoing rays which are asymptotically tangent to the 
same flow line of the Hamiltonian system induced by hf ^-\-V — E. The precise relationship 
between the integral kernel of Sh and the total sojourn relation is discussed in Section jS] 
and in particular we see that the canonical relation of Sh is a perturbation of the identity 
operator of order determined by a, the rate of vanishing of V at infinity. 
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We elaborate the latter remark in the special case that E = 1 and V is central {V{x) = 
F(|a;|)) on The bicharacteristics of the Hamiltonian p = + V — 1 are paths x{t) 

in R^ satisfying Newton’s equation x{t) = —2WV{x{t)). The crucial object related to the 
dynamical system in this context is the ‘scattering angle’ E, the angle by which an incoming 
ray is deflected by the potential [TH]- Indeed, in this case the canonical relation of is 
the graph of the map of T*S^ —> T*S^ taking a point (w, rj) to (w + E(? 7 ), 77 ), where (w, ij) 
corresponds to a straight ray Xo(t) = cot + r] and 77 T w. Here, the scattering angle is given 
explicitly by the formula [H Eqn. 2.6] 


E(77) =71 — 2 



_ V _ 

r2 


dr, 


where is the minimum distance to the origin of the bicharacteristic ray x{t) of the 
Hamiltonian p that is asymptotic to Xq for time near minus infinity. When V ^ c/r“, it is 
straightforward to compute that rm = 77(1 + 0{r]~°‘)) and E(77) = 0{rj~°^). 


This is the third paper in a series analyzing the asymptotic distribution of the phase 
shifts in the semiclassical limit using geometric microlocal techniques, the first two works of 
which consider smooth compactly supported potentials V [Si- It is instructive to compare 
the Main Theorem with the main result of [8] , which is 

Theorem 1.2. Let V be a real, smooth, compactly supported potential, and E gM. a non¬ 
trapping energy for the Schrodinger operator Elh = h'^A-\-V — E. Assume that the set of 
periodic points of powers of the reduced scattering map associated to Elh have measure zero 
in . Define the sequence of measures Vh on by 

(1.5) ^ 

Q^‘^f^rL,h ^spec(S/i) 

where I is the set of (01,77) £ associated to bicharacteristics that meet the support 

of V. Then, for every f £ supported away from 1, we have 

1 

(1.6) \imAvh,f) = —j /(e*"^)#. 

h->-0 271 7o 

In particular, the spectrum of Sh is asymptotically equidistributed on the unit circle away 
from the point 1. 

The main differences between the Main Theorem and Theorem 11.21 are 

• The rate of accumulation is different. There are about eigenvalues in a 

sector in the case of compact support [8] as opposed to for potentials 

decaying like 

• For compactly supported semiclassical potentials, the phase shifts equidistribute 
around the unit circle as h —^ 0, whereas for polynomial decay they do not; instead, 
we get the homogeneous distributions in the Main Theorem. 

• The equidistribution result for compactly supported potentials relies on two dy¬ 
namical assumptions on the bicharacteristic flow (the first being non-trapping). For 
polynomial decay, neither of these assumptions are required. 

These differences arise from the fact that Sh for a compactly supported potential is 
semiclassically equal to the identity operator outside a compact set in phase space, so the 
difference — Id is the difference of two semiclassical FIO’s of order 0 with compact 
microsupport. As such their traces grow like as h —>■ 0 (see e.g. O Appendix]), and 
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the volume of phase space on which — Id is microlocally nontrivial enters into the leading 
asymptotics. 

By contrast, in the present setting, compact subsets of phase space are irrelevant since 
their contribution to the measure /i?i is order x 

which is a positive power of h, so only the asymptotic behaviour of the dynamics is impor¬ 
tant. This explains why the nontrapping assumption is not relevant in the present setting, 
as trapped rays only occur in a compact region of phase space. 

The quicker rate of accumulation of eigenvalues, ^ as 0 can be 

understood heuristically by observing that, for each level of h, there is an ‘effective radius’ 
r{h) ^ outside of which the decaying potential V is 0{h), and therefore semiclassically 

negligible. This radius tends to infinity as a negative power of h, leading to an effective 
‘interacting’ volume of phase space that grows as . Since a unit of phase space 

volume contributes roughly eigenvalues, there are about 

phase shifts deflected away from 1. This observation also explains why we fail to have 
equidistribution, as one might naively guess based on Theorem 11.21 in the polynomially 
decaying case. Namely, there is no firm distinction between interacting and noninteracting 
parts of phase space, so correspondingly, there is no firm division between eigenvalues that 
are ‘essentially 1’ and ‘essentially different from 1’. So, unlike the compactly supported case 
where the measure divides into a finite, equidistributed part and an infinite atom at 1, in the 
polynomially decaying case, the point mass at 1 is ‘smeared’ into an absolutely continuous 
measure with infinite mass near 1. Thus equidistribution is not possible as it would only 
account for a finite amount of mass away from the point 1. 

One technical challenge of this work is that to treat potentials V for which one has only 
the derivative estimates of the main theorem requires an extension of the results in |12) . 
Indeed, the structure of the integral kernel of the scattering matrix, the Poisson operator, 
and indeed the outgoing and incoming resolvents are treated in |12] in the case that 

OO 

V = r-2^aj(a;)r-^ 

j=0 

for uniformly bounded aj G In other words, in the case that F is a smooth 

function of p = 1/r and w at p = 0. The potentials V under consideration here are merely 
conormal (see Appendix and thus some care is required to show that the scattering 
matrix has the FIO structure one would predict by analogy with the smooth case. This 
extension is done in detail in Appendix 

The introduction to [8] contains a literature review on the topic, to which we refer the 
reader. In particular, |20) contains an asymptotic formula for the phase shifts for central 
potentials of polynomial decay, in this case with two asymptotic parameters. Namely, they 
analyze /S.-\-gV for large g or k, in particular obtaining, when g = k'^ = l/h? a. formula 
for the phase shifts which implies our Main Theorem in the special case of central potentials. 
Other related work includes HHIHS]. 

2. The scattering matrix, Sh 

We now describe the FIO structure of the scattering matrix for potentials with polynomial 
decay. Note that by setting V = V/E and h = h/'/E we may reduce the general E case to 
the E = 1 case, and so we assume when convenient that 


E=l. 
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2.1. The canonical relation of the scattering matrix. We now describe the canonical 
relation of the scattering matrix. In fact, following [ 8 ], we define a Legendre submanifold 
of X X M related to the scattering matrix in a way we describe in detail in 

Section 12.21 

Given w' € and rj' G orthogonal to uj', there is a unique bicharacteristic ray 

of the seiniclassical Hamiltonian flow associated with Hh = 

H — 1 satisfying 

(2.1) = w't + T]' + o(l) for t « 0. 

and + V(xcj',n') = 1- (Recall that a (semiclassical) bicharacteristic ray 7 = (x,^) 

is a solution to Hamilton’s equations x{t) = 2,^, ^(t) = —\7V{x).) By our non-trapping 
assumption, as t —>• -|-oo, this ray escapes to infinity, taking the form 

( 2 . 2 ) Xui’,r]'{t) = Lo{t — t) +7] + 0 ( 1 ) for t » 0 . 

Here r = t{uj'^tj') is the ‘time delay’. As explained in [5], the pairs {uj',ri') and (^, 77 ) can 
be interpreted as points in The map (w'jjy') !->■ ( 01 , 77 ) is known as the reduced 

scattering map S = Se=i at energy E = 1. (For the arbitrary energy reduced scattering 
map see [S] of Appendix]^ The ray corresponding to (uj\r]') produces an additional piece 
of data, a function ip: —> R defined b}Q 

/ OO 

Xu:',r,'{s) ■ VV{Xu:’,r,'{s))ds. 

-00 

The reduced scattering map, together with the map (/?, determine a Legendre submanifold 
of X X R^, endowed with the contact form 77' • duj' + rj ■ dui — dip^ called the 

‘total sojourn relation’ in m- 

(2.4) L = {(o;', 77 ',o;,- 77 ,()!)) | ( 01 , 77 ) =S{uj’,r]’), </> = ip{uj’,r]’)}. 

That is, the contact form vanishes when restricted to L. This implies immediately that S 
is a symplectic map. The scattering matrix Sh can be described either as a semiclassical 
Lagrangian distribution with canonical relation given by the graph of S, or as a Lagrangian- 
Legendrian distribution associated to L, as we discuss further in Section [521 

We now find it convenient to switch to using local coordinates y = ( 1 / 1 ,..., yd-i) on the 
sphere We will use 77 for the dual coordinates on the cotangent bundle. This is a 

slight abuse of notation (compared to the usage of 77 above) but we find it convenient and 
should not cause confusion. For example, the contact form written in ( 7 /, 77 ) coordinates is 

11 ■ dy = J2iVidyi- 

Due to the decay of V at spatial infinity, the reduced scattering map S tends to the 
identity as | 77 '| —>■ 00 . It will be important in our analysis to understand precisely how 
this happens. In Appendix |A1 we will prove the following proposition which describes the 
behavior of the map S in the large 77 regime. Before we state the result, we recall the 
following standard terminology: we say that a function a = liy, 77 , h) is a symbol of order 
m in 77 , i.e. is in the space 5”’', if for each multi-index k £ Ng“^, k' £ Nq“^, there is Ck,k’ > 0 
such that 

(2.5) <Cfc,fc,(77)^-l'='l, 

where \k\= ki-\ -+ kd-i and {rf) = {\r]\‘^ + 1 )^/^ and Ck,k' is independent of h. 

^The function is closely related to, but not the same as, the time delay function r. See [S] Section 2] 
for further discussion. 
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Proposition 2.1. Suppose that a > 1, and that V S can he expressed in the form 

v="^ + w. 


where W satisfies m. Then the map {y,T],(j)) = {S{y',r]')) satisfies 


( 2 . 6 ) 


yt = y'i + a^{y',r]')\r|'\-°‘ + ei 
Vi = -V^ + bi{y\v')W\^~°' + e^ 
ip = ciy\p')W\^-" + e', 


for each i = — 1 and \ri'\ large, where ai,bi and c are smooth, Ci € S “ and 

ei,e' € 5'^“““'^ fsee (12.51 ) . ) 


2.2. The Schwartz Kernel of Sf,- Let us assume for the next few sections that if is a 
nontrapping energy, and defer the trapping case to Section [5] As shown in |12j . under this 
assumption, the scattering matrix is a ‘Legendrian-Lagrangian distribution’. This means 
that for each fixed h > 0, Sh is a (homogeneous) FIO; in our cas43, it is a pseudodifferential 
operator, in fact equal to the identity up to a pseudodifferential operator of order 1 — a. 
As h —>■ 0, the scattering matrix is, in each bounded region of phase space a semi- 

classical Lagrangian distribution, but it is not a semiclassical pseudodifferential operator. 
Instead, its canonical relation is the graph of the reduced scattering transformation, which is 
only asymptotically equal to the identity; Proposition 12. II makes precise how this happens. 
What is new about the result in [12] is that it gives the precise oscillatory integral form 
of Sh{E) in the transitional regime; that is, where the semiclassical frequency y tends to 
infinity, uniformly as /i —>■ 0. However, more regularity was assumed on the potential V in 
[12] : the assumption made there translates, in our context, to the potential having a Taylor 
series at infinity of the form explain how the result extends to the 

potentials considered here in Appendix lAl 

Remark 2.2. The reason that the term Legendre distribution is used in [12] is because Sh{E) 
is associated to the Legendre submanifold (12.41) . which gives an extra piece of information, 
namely ip, in addition to the symplectic map S. This eliminates an ambiguity (up to 
an additive constant) of the class of phase functions locally parametrizing the associated 
Lagrangian submanifold graph(5); see the appendix of [5] for more discussion on this point. 

We now express Sh{E) as an oscillatory integral microlocally near infinity. 


Lemma 2.3. Suppose E is a nontrapping energy. Then the scattering matrix Sh takes the 
form 

(2.7) Sh = Fi+F2, 

where Fi is a zeroth order FIO with compact microsupport, and F 2 has an oscillatory integral 
representation of the form 

(2.8) F 2 {y, y') = (^)"-' J + b{y, y', y, h)) dy, 


^The results of m apply to asymptotically conic nontrapping manifolds. In general, the absolute 
scattering matrix is a FIO associated to the canonical relation of geodesic flow (at inflnity) at time tt. In the 
case of one obtains the ‘relative scattering matrix’, which is what we consider here, by composing the 
absolute scattering matrix with the antipodal map and multiplying by This reduces the canonical 

relation to the identity in this special case. 
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where G{y',ri) is a symbol of order 1 — a in r], and b is a symbol in rj of order —a (see 
(ESI);. Moreover, 

( 2 . 9 ) G{y',r^) = g{y',f,M^-^ + ~g, g £ 


Proof. We first address the question of finding a phase function parametrizing L near infinity, 
that is, for large (77,77'). By LemmaEH we may use {y', 77) as coordinates on the Legendrian 
L in the large \r]\ region, and we write the remaining coordinates on L in terms of these as 

y = Wiy',r]), g' = N{y',g), ip = T{y',g), 


where {y' ,g' ,y,r], ip) £ L. The fact that L is Legendrian implies the following identities 
amongst these functions W, N, T, arising by expressing the vanishing of 77' • dy' + r] ■ dy — d(j) 
in these coordinates: 


(2.10) Ndy[ = 


d-1 

E 

7 = 1 


ar 


d-1 

E 

7=1 


17* 


dm 


, dT , 

drjj — ——drji = 0. 


dr]i 


Using these identities, one can check that the Legendrian L is parametrized by the function 

(2.11) $(7/', 7 /, 77 ) = iy-W{y',g)) • 77 + T( 7 /', 77 ). 

Let us write W{y',r]) = y' + W{y',rf). Then we have, comparing (12.9|) and (12.111) . 

(2.12) G( 7 /', 77) = -W{y', g)-g + T{y', 77). 

It thus suffices to note that by Lemma EH lU is a classical symbol of order —a, and T is a 
classical symbol of order 1 — a. Compare with [m Section 7.2]. 

That the scattering matrix has a local oscillatory integral expression using the phase 
function $ with symbol by a = 1 + 6 with b £ S~°‘ is shown in Appendix IA. 4 I □ 


Remark 2 . 4 . The choice to parametrize L using a function G{y' , 77) was an arbitrary one. 
We can just as well use {y,r]') to furnish coordinates on L, and then, writing 

7 /'= W'(7/,77') = 7/ +W'(y,77'), g = N’{y,r]'), ip = T\y,g’), 

the function 

^'(2/', V: V') = {y' -y) -r]' -W' -r]' + T'{y, 77') 

also parametrizes L. (To be clear, (p = T'{y, 77') means that T'{y, 77') is the value of (/? on L 
at the point (?/', 77', y, 77, ip).) Replacing the dummy variable 77' by 77 = —77', and setting 

G'{y,rj) = W'{y,-rj) - rj + T'{y,-Ti) 

gives = {y — y') -g + G'{y,Tf), and it follows as above that G' = g' {y' + g', g' £ 
gi-a-e ^ Furthermore, we claim that 

( 2 . 13 ) g'=g. 

In fact, since W{y', rf) = —W'{y, —rj) and T{y', 77) = T'{y, —rj), and writing rj = g +N{y',g), 
N £ we find that 

W{y',r]) = -W'{y' + W{y',g), -g + N{y',g)), 

from which it follows that 

^{y':V) = {y\ ~''l) modulo 

Similarly, T'{y',—g) — T{y',g) £ Thus we see that G — G' G from which 

( 12.1311 follows immediately. 
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Thus an oscillatory integral of the form in (12.81) with G as in Lemma 12.31 can also be 
written 

(2.14) + b{y, y', r,, h)) dy, 

where G' has the same properties as G, in fact G — G' £ and b has the same 

symbolic properties as a in (12.8|) . 

By adjusting the division between Fi and F 2 suitably, we may assume that G, as well as 
are sufficiently small, which we do without further comment. Indeed, choosing 
a function y = xi'd) with y = 1 1^1 < ^/2 and y = 0 for \r]\ > R and writing 

j ^^{v-v’)-ri+G{v',v)/h^-^ + a)dy = 

j + a) d77 + j - y)(l + a) dy, 

and taking R large enough and including the first term on the right in Tj produces the 
desired effect. We drop the y from the notation in F 2 since e.g. we could take a = —1 for 

\v\ < R- 


2 . 3 . Powers of the scattering matrix. To prove the Main Theorem, we will compute the 
trace of Sh{E)^ — Id for all integers k. Thus it is important to understand how to represent 
the powers as oscillatory integrals. 

First assume that fc > 1. The kth power of Sh is {Fi + ^ 2 )^, and if we expand this 
product, every term has compact microsupport except for F^. (See Section |4] for a further 
discussion of the other terms in the expansion and why their contribution to the eigenvalue 
asymptotics is lower order.) For fc = — 1, recall that = 5^, and thus the integral kernel 
of is given by the hermitian conjugates Fj* +F| in ()2.7|) . Here F 2 (y, y') = F 2 {y', y). By 
Remark |2.41 we may take the phase function of F 2 to be of the form in ()2.I4|1 . specifically 
with G' as in (I2.I4|1 . we have 

(2.15) F*( 2 /, y') = F 2 {y', y) = (^)"-' J + b{y,y', -y, h)) dy, 

We will show that F 2 has the following oscillatory integral structure. 


Lemma 2.5. Suppose fc > 1. Then the FIO F 2 has an oscillatory integral representation 
of the form 


(2.16) 


< 2 ^> 


d-1 


gi{{v-v'yv+kG{y’.v)+E^{v'.ri))/h^^ ^ 


where Ek is a symbol of order 2(1 — a) in y, and bk is a symbol of order 1 — a. Similarly, 
(E^)^ has an oscillatory integral representation of the form 


(2.17) 


+ b,{y, y', y, h)) dy, 


Proof. See Appendix IB] 


□ 


The key point here is that, up to a term vanishing faster as \y\ —>• 00, the effect on the 
phase function $ of raising the scattering matrix to the power fc is essentially to replace G 
by fcG. 
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3. Proof of the main theorem 


The main idea of the proof, motivated by [131 HI] and [5] , is the following observation: 
if {vh} is a family of finite measures on parametrized by h > 0, and if each Fourier 
coefficient of converges to that of a certain finite measure u as /i —^ 0, then converges 
to u in the weak-* topology. In our case, however, we cannot apply this directly, as the /x/j 
in (11.31) are infinite measures. Instead we have the following variant. Consider the following 
weighted sup norm for functions on 

(3-1) ll/IL = sup 

and the associated Banach space 

(3.2) = {/ e :3gG such that f{z) = (z - l)g(z)}. 

Then we show 


f(.z) 

z-1 


Proposition 3.1. Suppose that the (infinite) measures Vh o,nd v act (by integration) as 
bounded linear functionals on C')),(§^). Moreover, assume that the norms of Vh in the dual 
space are uniformly bounded in h. Then if 

(3.3) lim f p{e'^^)dvh = [ p(e''^)dv 

^“^0 dsi Jsy 

for all polynomials p € C')),(S^), then Vh ^ v in the weak-* topology on every compact subset 
o/§i\{l}. 


Proof. The proof of this proposition consists of two elementary steps. We first observe, as 
in [31 Proof of Lemma 5.3], that polynomials in (7° (S^) are dense in C° (S^). The proof is 
so simple that we repeat it here: given / G (7° (S^), by definition f = [z — l)g(z) for some 
g G (7°(§^). We approximate g in (7°(S^) by a sequence of polynomials pj. Then {z — l)pj 
lie in C'°(S^) and approximate / in the (7° (§^) norm. 

Then, given e > 0, and a continuous function / on the circle, supported away from 1, we 
need to show that 



< (7e, 


provided h is sufficiently small. We choose a polynomialp in (7° (S^) such that ||p—/||[,„] < e. 
Then we estimate 


(3.4) 





fdvh 


The first term is bounded by Ci\\p — /||[u,] where (7i = ||u||(co). is the dual norm of v. 
The second term is bounded by e provided that h is sufficiently small, using (13.3L The 
third term is bounded by C2\\p — f\\[w] where C2 is a uniform bound on ||u;j||(c' 0 ).. Taking 
(7 = Cl + (72 + 1, this completes the proof. □ 


In the case of interest, Vh will be the measure gh dehned in (II.3|) and v will be the 
pushforward of a homogeneous measure. In view of Proposition 13.11 we need the following 
result. 


Proposition 3.2. There exists c > 0 such that for h sufficiently small 

(3-5) \{l^h,f)\<c\\f\\^. 
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We also need 


Proposition 3.3. For every fc G Z, — Id is trace class. 

Propositions 13.21 and 13.31 will be proved in Section [5l 

It is easy to see that the polynomials — 1, for fc 7 ^ 0 G Z, form a basis of the polynomials 
that vanish at 1. We thus need to show that, for each k, the quantity 

(3.6) J - l)d^fc = Trace(S'^ - Id) 

converges as h —>■ 0, and to find a limit measure /r such that the limit of (13.61) is equal to 


(3.7) 



l)dy.. 


It turns out that the limit of the quantity (13.61) is given by a power c±|A:|'’', where the 
coefficient c± depends only on the sign of k. Such a homogenous ‘Fourier series’ comes from 
a ‘homogeneous measure’. We now describe precisely what this entails. 


Definition 3.4. Let j3 < —I. We say that a measure fi on \ {1} is the pushforward of a 
homogeneous measure of degree fi if there is a measure on R 0 / the form 


ci0Pd6, 0 > 

C2\0\^d0, 0<O 


such that fi is the pushforward of v under the quotient map P : R 1 —>■ = R/27rZ.' 

p, = F,(i/) on \ {!}. 


To state the following lemma, we will need to define the constant 


(3.9) 



l)0-^-^d0. 


for 0 < 7 < I. Note that F is finite for 7 in this range. 


Lemma 3.5. Suppose that p is a measure on that lies in the dual space of and 

is such that, for some 7 > 0 , 


(3.10) 



(Fci+rc 2 )fcT', fe = 1,2,3,... 

(rci+rc 2 )|A:|T', fc = -I,-2,-3,... 


Then on \ {!}, p is the pushforward of a homogeneous measure of degree f = —I — 7 . 
Indeed, it is given by v in ()3.8I) (with the same constants ci and 02 ). 


Remark 3.6. In this lemma, we are only making a statement about p away from the point 
1. Notice that there could be an atom at 1 about which we can say nothing, as this would 
not affect the integrals in (13.101) . 

Proof. We first note that the integrals in (13.101) determine p uniquely as an element of 
the dual space of C® (S^), hence uniquely as a measure away from the point 1. This is an 
immediate consequence of the density of polynomials in C')),(S^). In view of this, it suffices to 
show that the measures in Definition [331 homogeneous of degree fi = —1 — 7 , have ‘Fourier 
coefficients’ of the form ([3.101) . 
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With /3 = —7 — 1, let V = H{9)6^d9 where H is the Heaviside function, and let ^ be the 
pushforward of v. We compute, for fc > 0, 


Jsi Jr 


(3.11) 


- 1)9^ d9 


= /c^ / (e*®- 


= k^T. 


Similarly, for A: < 0 we Hnd 

[ -l)d^l=\kp [ {e-^^-l)9>^d9 = \k\'^T. 

Jsi Jo 

Similarly, if u = (1 — H{9))\9\^, then for fc > 0 

p pO poo 

/ (e*'=^-l)d/i= / -l)\9fd9 = / - 

Js^ J — 00 J 0 

and /gi — l)d^ = [fc^F for fc < 0 . 

We can now prove the Main Theorem. 


= fc'^r, 


□ 


Proof of the Main Theorem. Proposition 13.21 shows that the measures fj.h are uniformly 
bounded in the dual space of (7° (§^) (since a > d). This means that we can apply Proposi¬ 
tion [ST] showing that jih converges to a measure /i on provided that the convergence on 
polynomials vanishing at 1 in (13.31) holds. Since the polynomials — 1 for fc ^ 0 G Z are a 
basis for polynomials vanishing at 1, it suffices to check ()3.3|) for these polynomials. So this 
requires computing the limit, as fc —)► 0, of Trace(S'^ — Id). This we shall do in the Section 
m and the result is (|T6l) and (ITSI) . That is, the integrals are given by ^ckP 

for fc > 0 and {yfE/2TT)'^~^ck^ for fc < 0, where c is the constant in (14.111) . in particular 
c = oiP -I- a 2 r for 01,02 in (I4.11|) . Thus by Lemma TS.51 the homogeneous measure ^ given 
by the pushforward of ai{'/E/2 ttY~^H 9^ a 2 {'/E/2t:Y~^{1 — H)\9\^ which pairs with 
gifes _ gjyg above values. □ 

We now prove Corollarv ll.il 

Proof of Corollary \1.1[ Given 0 < (j)o < (j)i < 2 tt, and let l[ 0 o, 0 i] • ^ ® Ih® indicator 

function of the corresponding sector of the circle, l[ 0 o, 0 p(e®®) = 1 if (l)o < 9 < (j>i modulo 
271 and is zero otherwise. Then for N{(j)Q, (fi) as defined in the corollary 

= Tr(l[0g_0p(S';t)) 

Let / and g are continuous, non-negative functions, on the circle supported on \ 1 such 
that / < 1 ( 00 . 0 i] < g. Then 

Tr/(5i.) < Tr(l[0o,0o](^/.)) < T^9{Sh), 
and all three quantities are finite. Thus 

y /fc + 0 ( 1 ) = (fc/i, /) < h°‘'>N{(l)o, (fi) < {gh, 9 ) = J 99 + o(l), 
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where o(l) denotes a quantity which goes to 0 as h 0. But for any e > 0 we can choose / 
and g can be chosen so that J fg. = J — e, J gg = f + e, and thus for any 

e we have 

J ^l(j>o,(j>i](^ ~ ^ + o{l) < < J l[ 0 jj + e + o(l), 

so limh“'>' 7 V(^ 0 j </>!) = / l[ 0 o, 0 i]^; proving the Corollary. □ 


4. Traces and compositions 


We now compute the traces 

(4.1) lim/i“'^Tr(5^-Id), 

h—^0 

again with 7 = (d — 1 )/(q; — 1), assuming still that a > d. Indeed, we will prove 


Lemma 4.1. There is a constant c such that for fc € Z, 


lim(^^,z'= 

h —>^0 


J c( 27 r)-(‘^-i)fcT' 

\ c(27r)-(‘^-i)|fcp 


if k >0 
if k < 0 


Indeed, c = aiT + a 2 r where oi and 02 are defined in (j4.11|) 


Note that any semiclassical, zeroth order FIO with compact microsupport has trace 
bounded by (see e.g. [5J Appendix].) Therefore, in the limit above we may replace 

S’f by F^. Indeed, in 

fc-i 

Sf-ld = F^-ld + Y^ 

j=o 



all the terms in the sum on the right have compact microsupport, are thus trace class with 
trace bounded by and as we will see, the trace of Ff — Id increases at the rate 

Since 


= a{d — 1)/{a — 1) > d — 1, 

the trace of F 2 — Id contributes the leading order part of Tr(S'^ — Id). For the same reason, 
we can replace the identity operator by another FIO that differs from it by an operator with 
compact microsupport. So we can restrict to the microlocal region ji^j > R for arbitrary 
R, and using = h'^ and Lemma 12.51 we can write the terms Schwartz kernel of 

— Id which contribute to the trace to leading order in the form 


(4.2) 


J + dr?- J dg, 

\ri\>R hl>R 


We then compute ()4.1I) by setting y = y' and integrating over y. To be more precise, 
the Schwartz kernel corresponding to the oscillatory integral expression in (14.211 is actually 
a half-density on x acting on half-densities on meaning that if, for the 

moment, H{y,y') denotes the distribution in (14.2|) . then id(i/, 2 /')|d 2 /dy'|^/^ is the Schwartz 
kernel - more precisely the Schwartz kernel is a finite sum of these - and it acts on half 
densities (j){y)\dy\^^'^ by 


( 4 . 3 ) 


^ ( / H{y,y')(l){y')\dy'\)\dy\^^^. 




DISTRIBUTION OF PHASE SHIFTS FOR POTENTIALS WITH POLYNOMIAL DECAY 


13 


It is standard that the trace of this operator is / H{y,y)\dy\, and thus we are tasked with 
computing 


(4.4) I e^>^^^y’^^/^akiy.y,v,h)dydy+ j - l) dydy 


lv\>R 


\ri\>R 


Since at = 0(|i?|^~“) and a > d, the first integral is absolutely convergent. Due to the 
positive power of h out the front, this term is zero in the limit /i —>■ 0 . 

So consider the second term. We write 

(4-5) G{y, y) = \y\^~°‘g{y, y) + \y\^~°‘~^Ky^ v), 

where 5 is a symbol of order 0. We change variable to y' = y{h/kY/^°‘~'^'> to obtain 

h~i 


{2ttY 


J (^e^i^G(v,v)/h 

R 


\v\>R 

fcT- 

( 27 r) 


The integrand in this integral is dominated for small h by C\y'Y~°‘ for C > 0 independent 
of h, which is integrable as a > d. Thus, by the dominated convergence theorem, we can 
take the pointwise limit inside the integral, and obtain 

(4.6) lim(^,„z'=-l) = lim/i(‘'-i)“/(“-i)Tr(,5^-Id) =c(l/(27r))'^-iF, 

h^O /i—>0 


= f dydy', k>l. 


where 

(4.7) 

The fact that Sh is unitary immediately implies 

(4.8) lim-Id) = c(l/(27r))'^-i|fc|'>' fc <-1. 

h—^0 

It remains only to evaluate c. Write 

9{y, v') = 9+{y, v') - 9-{y, v'), 

where g-^. = max{( 7 , 0}, = max{—g,0}. Then 

J ^ J - 1 ) dydy' 

+ f - 1) dydr/. 
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Considering the first term on the right hand side, we write rj' = rr]' with 77 ' = 77V 1 ^^ I 
for fixed y,fi' with g+[y,f}') ^ 0 , compute 



where in the first step we set r = in the third step we set “ = p, and where F 

is the constant defined in (13.91) . A similar computation shows that 


(4.10) 



-ig-{v,fi’)W\ 




gldfj'dy. 


and thus it follows that 

(4.11) c = oiF + a 2 r, where oi 



g'J_dfi'dy, and 02 



gldfj'dy 


5. The asymptotic distribution of phase shifts 

The aim of this section is to prove Propositions 13.21 and 13.31 to which end we must first 
obtain an estimate for the rate at which eigenvalues of Sh accumulate at 1 

5.1. Eigenvalue accumulation at 1. 

Proposition 5.1. Let be the eigenvalues of the scattering matrix Sh- There 

exists a constant c > 0 such that for each e > 0 sufficiently small, 

(5.1) #{n : - 1| > e} < ce-''h-°'^, 

where as above 7 = (d— 1 )/(q:— 1 ). 

Before we prove the proposition, we remind the reader of the following fact; let Kh be a 
semiclassical pseudodifferential operator of semiclassical order 0 on a compact manifold N 
without boundary, with compact microsupport. Then there exists a c > 0 such that off a 
subspace Wh C Lf with diurWh < we have 

(5.2) WKhWwf^L^ = 0{hn- 

Indeed, this can be shown using properties of the semiclassical Laplacian hf/S-M correspond¬ 
ing to a Riemannian metric on N as follows. If one considers a smooth, compactly supported 
function y: T*N —>• R, with WYh{Kh) C suppy, then writing 

Kh = Khx{h^K) + Kh{ld -y(/i^A)), 

the term on the left hand side satisfies ||AAy(h^A)|| < ||Ar 7 i||||y(h^A)|| where the norms 
are operator norms as maps on Lf. (see [25] for a definiton of the semiclassical wavefront 
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set WFh of a semiclassical FIO.) By the Weyl asymptotic formula for semiclassical pseu¬ 
dodifferential operators of order 2 [25l Section 6.4], which states that, if are the 

eigenvalues of A 

(5.3) #{j : Af < A^} = -h 0(A'^™^-^), 

as A —>■ oo, we see that x(h^A) is identically zero off a set of dimension no bigger than 
cVol(supp(x))/i“‘^“^^. For the term on the left. Id—x(h^A) is a semiclassical pseudo of 
order 0 , with microsupport in (suppy)"^, and thus 

WFhiKh{ld-xih^^))) C WFh{Kh)Fi (suppx)‘= = 0, 

and in particular ||iF/i(Id —x(h^A))|| = 0{h°°). 

In the proof we will use a semiclassical version of the Calderon-Vaillancourt theorem, 
which, in the non-semiclassical setting [TJ Theroem 2.73] states that a pseudodifferential 
operator Q on R", 

Q = J e(^-^')<a(z,z',C)dC, 

satisfying 

||a|| 2 n+i := sup ||5“^,a^a|]L<x, < oo 
|a| + |/3|<2n+l 

is bounded on with ||( 5 ||L 2 _,,i 2 < c„|]a|| 2 n-i-i, where c„ is a constant depending only on 
the dimension n. Setting ( = C/h shows that for a semiclassical pseudodifferential operator 
of order 0 , 

Qh = h-^J 
if 

(5.4) ||5||2n+i.ii := sup \\d^^,{hd(;)^b\\L^ < oo 

|a|-|-|/3|<2n-|-l 

is bounded on with norm bounded by c„|| 6 || 2 n+i,/t. Indeed, this is just Courant-Vaillancourt 
with symbol depending on a smooth parameter applied to a semiclassical symbol 6 (z, z\ h) 

Proof of Proposition I5.ll Consider the operator Sh = F 1 +F 2 decomposed as in ( 1 ^ . where 
Fj has compact microsupport and F 2 consists of a finite sum of terms of the form ( 12 . 8 F 
We begin by taking a cutoff function xr- —>■ M with xr{v) = 1 for \r]\ < R/2 and 

suppx C {|? 7 | < R}. Taking R sufficiently large we have 

F/, - Id = {Sh - Id) Op;,(xr) + {Sh - Id)(Id - Op^(xr)), 

where Oph{Q) applied to a symbol q G 5'™(S‘^“^) denotes the right quantization of q to 
a semiclassical pseudodifferential operator of order m, again see [25] . Here, the operator 
{Sh — Id)Op;j(x) G and thus the statments proceeding the proof apply, and since 

OP/i(xr) is a semiclassical pseudodifferential operator of order 0 with compact microsupport 
[U Appendix], again using d — 1 < ay, we see by the discussion proceeding the proof that 
the behavior of {Sh — Id) Op,j(x) has no bearing on (15.11) . On the other hand, we can take R 
large enough so that WF/i(Fi Op/j(xr)) = 0 and thus the Schwartz kernel of the operator 

Ah = {Sh - Id)(Id - Op;,(xr)). 

is a sum of terms of the form (12.81) plus terms of order 0{h°°), and we thus focus our 
attention on terms as in ( 12 . 8 L 
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By Remark 12.41 we may take F 2 to be as in (I2.14|) with phase function $ = {y — y') ■ 
T] + G{y,ri) which appear in F 2 have the asymptotics G = a{y,rj)\r]\^~°' + 0(|77|““). For 
constant ^ > 0, we consider two asymptotic regimes 

regime I: > S here e:}q)iG{y,r])/h is not oscillatory 

(^■^) i/n- -1 

regime II: \r]\h '^ “•' < 5 here expiG(y,7y)//i is oscillatory. 

As in the discussion proceeding the proof, we use functions of the semiclassical Laplacian 

Ph ■= h^Agd-i. 

Let y: R'*’ —R be a bump function with y(r) = 1 for r < 1/2, y > 0, and suppy C [0,1], 
and write 

( 5 . 6 ) Ah = Ah^i + Ah.2i 
where 

Ah,i ■■= A^y((e/i)i/(“-i)p^), Ah,2 := +A,,(Id -y((eh)i/(“-i)p^)). 

We analyze these two operators separately. 

For Ah,i^ we begin by pointing out that there exist subspaces I 4 C with dimR < 
g- 7 ^-a 7 ^ Indeed, by the Weyl asymptotic formula 

(El, the direct sum of the eigenspaces of ^'>Ph with eigenvalue less than 1, which 

we take to be I 4 , satisfies 

(5.7) dim 14 < = ce-^h-°^^. 

Thus Ah,I is also identically zero off Vh- 

Now consider Ah ,2 ■= A^,(Id —y((eh)^/(““^)P;,)). We claim that there exists a constant 
G > 0 independent of h and e such that 

(5.8) ||A^,2|| < Ge. 

The Schwartz kernel of the operator y((eh)^/(““^^P;i) is given by finite sums of terms 

(5.9) (27r/i)-(‘^-i) j 


where y(y, 2 //? 7 ,h) is a semiclassical symbol of order zero with y||^|>i = 0{h°°). The 
operator Ah ,2 is given by terms of the form 

(27rh)-(''-l) [ F(V-V’yFh + iGiy,^)/h^^y^ y,^ _ ~^y^ ( £ ^ ) ^ Z ^ “ " ^ ) , ^ ) ) di? 

(5.10) ■' 

- (27r/i)-(‘'-i) / p(^-^')-''/"(l-y(y,y',,7(eh)i/(“-i),h))dr;. 


This requires some explanation. To compute the composition we must compose an operator 
whose Schwartz kernel is an oscillatory integral as in (I2.14|) . call it //i(y, y') with an operator 
whose Schwartz kernel is an oscillatory integral of the form (15.9|) . This is done by arguing 
along the lines in Appendix [B1 where in particular we see that the composition of two such 
oscillatory integrals is given by / Ih{y, y")Ih{y'', y')\dy''\- The situation here is substantially 
simpler since the operator on the right is a semiclassical pseudo, and the expression above 
is obtained easily from stationary phase. (Again, see Appendix IbI') 

Setting rj = h = (eh)^Z(“-i)/j factoring gives 

J ^■^Z^d(y, y', rj, h)drj, 


(5.11) 






DISTRIBUTION OF PHASE SHIFTS FOR POTENTIALS WITH POLYNOMIAL DECAY 


17 


where, letting h remain as a function of h and e for the moment, 

a{y,y\v,h) = (exp{iG{y,7j{eh)~^/^°~^^)/h)b{y,y’ 

X (1 - x{y,y',rj,h)). 

The immediate effect of this change of variables is that 

(5.12) F = exp(iG(i/, x (1 - x{y, v', V, h)) 


satisfies F — {1 — x) S “, i.e. it is a semiclassical symbol of order 1 — a as a function of 
rj and h (or rj and h for that matter.) Indeed, recalling that support of 1 — y is contained in 
{rj > 1} and thus |? 7|“^(1 — x) is bounded below for any 5 > 0, we claim first that, notation 
as in (|4.5L 


{^-x)G/h={l-x)i9{y,r])e\v\^ “ + 5 ( 1 /, y', 77(e/i) ^^)/h) 

satisfies that (1 — x)y//i G as a function of rj. Indeed, 

(5.13) I3,va|(i-x)9/'.|<ci.-'l 1 : a.V'C’h-aa.va??!, 


and while for /3' ^ 0,d~ ^ (1 — x) is compactly supported in y, the symbol estimates for g 
give that for /3 0 

(5.14) 

1(1 - X)dl9\ < (1 - 

< (1 - ^)(e/i)i/(i-«)-(i-«-=-l/3|)/(«-i) 

where in the last line we used that 1 — y is supported in \rj\ > 1, while for /3 = 0 firstly that 

1(1 - my,y', {eh)-^/^^-^'>/h\ < C{1 - x)(i?(e/i)-'/^“-'V-“-^ < Ce. 

The estimates in (15.131) and (15.141) together show (1 — x)Glh and thus f — (1 — y) € eS^~°‘. 
Moreover, 

a(y, y', h) = 1 + b(y, y', h) 

remains a symbol of order 0 and b a symbol of order 1 — a, and as in the case of F, 
b(y,y',rj{eh)~^/^°‘~^\h) G inly. Thus 

a=iF-il-x))+Fb 

is e times a semiclassical symbol of order 1 — a whose derivatives in y, y' and rj are uniformly 
bounded. In particular Calderon-Vaillancourt (see (15.41) and below) gives (15.81) . 

To hnish the proof, given e we take t' = ejC with C in (j5.8p and use e' in the arguments 
above to see that in Ah = Ahx + Ah .2 (see (15.61) 1. Ah.i has at most eigenvalues at 

distance e from 1 while Ah .2 is norm bounded by e. This is exactly the desired result. □ 
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5.2. Proof of Propositions 13.21 and 13.31 

Proof of Proposition \S.‘A Let p gN. Let 

Ah{p) = € specSh : 2-^p-^^ > - 1 | > 2 -^}, 

where elements are included with multiplicity. Then taking e = 2~p in (15.11) gives 

(5.15) \Ahip)\ < 

The pairing of / with is given by 

OO 

(5.16) W, /) = E /(e"*'"'*-") = E( E /(e"*'"'*-"))- 

specSh P=0 Ah{p) 

But, 

I E fi^^'M\<\\f\L E le"*'"'*’"-!! 

Ahip) Ah{p) 

(5.17) <\\f\L2-P\A,{p)\ 

< c ll/IL 2-P(2P("*-^)/(“-^)/i-“/(“-i)) 

< c ll/ll 

Thus 

(5.18) 

OO OO 

{ph, f) < Il/||„ ^ < c ||/| 1 „ E 

P—0 P—0 

and the above is summable if and only if 

1 ^-1 

1 > -- a > d, 

a — 1 

which is exactly our assumption on a. □ 

Proof of Proposition \3.A It is enough to prove for fc = 1, as for any other value of fc, we can 
write — Id as the product of Sh — Id with a bounded operator. For /c = 1, Proposition [ST] 
shows that the number of eigenvalues z (counted with multiplicity) of Sh such that 1 ^ — 1 | G 
[2~\2~^~^^] is bounded by Ch~°‘'^2^'^. Since 7 < 1, we can sum 2“'1+^ x Ch~°‘'^2^'^ over 
j G N, which is a bound for the sum of |z — I| over all eigenvalues z. It follows that the 
trace norm of Sh — Id is finite. □ 


6. Trapping energies 

Suppose now that i? is a trapping energy for the potential V. In this case, we write the 
scattering matrix ShiE) as the scattering matrix Sh{E) for a different potential V, which 
is nontrapping at energy E, plus a small remainder. We can choose the potential V to be 
equal to V near infinity. To do this, we first choose a function (f G C'“(K.+ ), equal to 1 in 
a neighbourhood of 0, and monotone nonincreasing. Then V := V + 2E(l){\x\/ R) will be 
nontrapping at energy E, for sufficiently large R. 

We then express the scattering matrix Sh{E) in terms of Sh{E). To do this, we follow 
[U Section 8B]. Let Rh = {h‘^A + V — {E + i0))~^ and Rh = (h^A + t^— (£’ + i0))“^ be the 
outgoing resolvents for the unperturbed and perturbed potential, respectively. Also, let Xiy 
i = 1,2,3 be cutoff functions supported near infinity in R", equal to 1 for |a;| > 2R and 0 
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for \x\ < R, such that XiXj = Xj when j < i. Then, following the derivation of [HI Equation 
(8-7)], i.e. taking H = A -|- E and A = Ejh? in that equation, we obtain 

( 6 . 1 ) X2RhXi = X2RhXi + X2Rh[x^, V]Rh[h^^ + V,X2]RhXi- 

These Schwartz kernels are defined on x x (0,/io]/t. As discussed in [12], if we use 
polar coordinates x = (t, w), x' = (r',w'), multiply the kernel of Rh by and take 

the limit r' —>■ oo, we obtain the Poisson kernel Ph{E), which is a function of {x,lo', h). If 
we then multiply the Poisson kernel PhiE) by and take the (distributional) limit 

r —^ oo, we obtain the kernel of the absolute scattering matrix; multiplying by and 

composing with the antipodal map A, we obtain the scattering matrix Sh{E) as we have 
normalized it. The same operations applied to Rh produce Ph{E) and Sh{E). Applying 
these operations to mi, we obtain 

(6.2) Sh{E) = ShiE) + i^^-^^/^AP;:{E)[x3,h^A + V]Rh[h^A + V,X2]Ph{E). 

For brevity, we write this in the form 

(6.3) ShiE) = ShiE) + BhiE); 

clearly Bh{E) is a uniformly bounded family of operators on L^(S‘^“^). 

Our previous arguments apply to Sh{E), since if is a nontrapping energy for V. So it 
suffices to show that the perturbation Bh{E) has no effect on the weak-* limit fl of the 
measures jxh associated to Sh{E), as h —0. 

To show this, we now cut off to small and large frequencies using a cutoff x(/i^Asd-i), 
where x(l) is compactly supported, and identically 1 near t = 0. For simplicity we write 
this operator simply as x- Thus we write 

(6.4) Sh{E) = xShiE) + xBh{E) + (Id -x)Sh{.E) + (Id -x)Bh{E). 

The first term is an FIO with compact microsupport, hence has trace norm bounded by 
The second term also has trace norm bounded by since this is true 

of X which is also an FIO with compact microsupport. The third term is the principal 
term, and the fourth we bound using wavefront set results. In fact, according to m, the 
semiclassical wavefront set of Ph{E) is contained in 

{(w, Tj; X, I the bicharacteristic through {x, has asymptotic t^ri + u}(t — to),t ^ oo} 

when the point x is restricted to a fixed compact set. Now consider the composition 
{ld—x)APfl{E)[x 3 ,h^A + V]. Composition on the right with [x 3 ,h^A + V] restricts the 
wavefront set to points x € suppVxa, that is, to x lying in some fixed compact set in 
R"^. On the other hand, composition on the left with (Id—x) restricts the wavefront set to 
points (oj, rf) in the support of the symbol of %. By choosing x suitably, we can arrange that 
this support is contained in \r]\ > R' for R' arbitrary. By choosing R' sufficiently large, we 
arrange that the wavefront set of (Id —x)AP^{E)[x 3 , h^A -|- V] vanishes. That implies that 
the Schwartz kernel of this operator is smooth and 0{h°°). The trace norm of the fourth 
term in (I6A1) is therefore 0{h°°). 

Now consider all the terms in Sh{Ep — Id, where Sh{E) is decomposed according to ()6.4() . 
The main term, ((Id —x)*S/t(A))^ — Id, is treated as in Sections El-El All other terms have 
trace norm bounded by and therefore their contribution to h'^°^ Trace(/S'?i(i?)^ — 

Id) vanishes in the limit h —0. This completes the proof of the Main Theorem in the case 
of a trapping energy. 
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Appendix A. Regularity of the sojourn map 

In this appendix we prove Proposition 12.II and Lemnia l2.3l Our first task is to determine 
the regularity of the Legendre submanifold L (12.4|) as |7?| —oo. To do this, we use the fact 
that L is the boundary value of a Legendre submanifold SR over a space of dimension one 
greater which is a bicharacteristic flowout, that is, the union of bicharacteristic rays. We 
start by defining some spaces of conormal functions, and then proceed to describe SR and 
its ambient contact manifold. 

This process will use the language, developed by Melrose [lSJ[T6], of analysis on manifolds 
with corners. Though some of this is quite involved we will provide some brief explanations 
and definitions for the convenience of the reader. 


A.l. Conormal regularity of solutions to ODEs. Let M be a manifold with corners, 
with boundary hypersurfaces Hi,... ,Hm and boundary defining functions pi,. ■ ■ ,Pm re¬ 
spectively [TS]. Thus the boundary dM is equal to the union of the Hi, and for each i. 
Pi is a non-negative, smooth function on M with Hi = {pi = 0} and dpi ^ 0 on Let 
p = Pi..., Pm be the product of boundary defining functions. We say that a vector field V 
on M is a b-vector field if it is smooth, and tangent to each boundary hypersurface Hi, or 
equivalently if V{pi) = 0{pi) for each i. 

Let e = (ei,..., em) £ R"* be a multiweight, one for each boundary hypersurface of M. 
The (L'^-based) space of conormal functions with weight e, A-{M), is defined as follows: 


(A.l) AHM) 


f € C°°{M°) 


^-f e L°°(M) and p-^Vi ... Vfc/ £ L“(M) 
for any k b-vector fields Vi,..., Vfe 


Here p^ is shorthand notation for the product p\^ ...p^. That is, / £ p^L°°{M), and 
remains in this space under repeated differentiations by b-vector fields on M. The space A- 
is a Frechet space whose metric we describe below for a simple example. 

We also use the notation 


= C°°{M) -f A-{M). 

The regularity condition of our potential V in the Main Theorem can be phrased in terms of 
the above spaces; the assumption on V can be expressed in terms of the radial compactifica- 
tion of R.'^, where 1/r is taken as the boundary defining function at the ‘sphere at infinity,’ 
and is equivalent to assuming that for some 0 < e < 1, £ C°°’*‘(R‘^). (Equivalently, 

V £ r““C'°°’'‘(R‘‘*).) We abuse notation slightly by defining a smooth map u = {ui,..., u„) 
on M° with values in C" to lie in A-{M) if and only if its components do. 

We note, for later use, the following result. The proof is straightforward and omitted. 

Lemma A.l. 

(i) If f £ is bounded away from zero, then 1/f £ 

(ii) If S : M ^ N is a between manifolds with corners M and N such that all 

components of S have regularity C°°'-{M), and S is invertible in the sense that it is invertible 
as a map and its Jacobian determinant is bounded away from zero, then the inverse map 
has regularity C°°'-{N). 

(Hi) Let 7i,...,7m be positive exponents, and suppose that e > 0 is sufficiently small 
(relative to the ^i). Then the statements (i) and (ii) above also hold if the space C^’- is 
replaced by C°° + Ilp]'C°°’-. 


^This means that the inverse image of every boundary defining function on Al is a product of boundary 
defining functions on M, times a smooth non-vanishing function. An invertible b-map induces, in particular, 
a bijection between the codimension fc-faces of M and the codimension fc-faces of N. 
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It is well known that solutions of ODEs 

^ = 2 /( 0 ) = 2/0 

are smooth if F is smooth, and y also depends smoothly on the initial condition yo- See e.g. 
Hartman m Chapter 5 ]. Here, we note the following variant of this standard result. We 
find it convenient to write the ODE in terms of a b-derivative, tdt- 


Proposition A.2. Consider the ODE 

dz 

t— = F{z,s,t) 2(0) = zo, 

(A.2) f 

' ' ds 

t— = sG{z,s,t) s(0) = So > 0. 

for 2 G and s G K+. 

(i) Suppose that F,G € x 1R+ x K^), Pi > 0, where /3i refers to the s variable 

and P2 to t. Then the solution 2 = 2(2o,so,t),s = s(zo,so,t) satisfies 

2(20, SQ,t) - 2 o, —(s( 2 o,so, 2 ) - So) G X M+ x Rf) 

sq 

locally near 2 = 0. 

(a) Suppose that F,G € x R+ x R^) + (Rf x R+ x R^), Pi > 0. Then the 

solution 2 = 2(20, So, 2), s = s( 2 o, so, 2) satisfies 

2 - 20, -(s - so) G 2C'-(RP^ X R+ X R+) + x R+ x R+) 

^0 

locally near t = 0 . 

(Hi) Let Pi = ^i + e, where 7^ > 0 and e is sufficiently small. Suppose that F,G ^ 
20 °° (RP X R+ X R+) + sT'D'^2C'°°(RP x R+ x R+)+A^ 1 ’^=(Rp x R+ x R+). Then the solution 
z = 2(20, So, 2), s = s(2o. So, 2) satisfies 

2-20, -(s-so) G 2C°°(RP„ xR+ xR+) + s2D^=C“(RP„ xR+ xR+)+A''^’^nRP„ xR+ xR+) 
locally near 2 = 0. 


Proof. We start by making some reductions. We first let 2 = 2 — 20 and s = log(s/sQ). 
Then 2(0) and s( 0 ) solve the initial value problem 


(A.3) 


dz , 

2 — =F(2 + 2o,soe , 2 ) 

0 

II 

0 

ds 

2— =G(2 + 2o,soe®,2) 

s( 0 ) = 0 . 

J) into a new variable Z, satisfying an 

dZ 


2 — = H{Z,zo,so,t) 

0 

II 

N 


and show conormal regularity in the (so,2) variables. 

Let S', T denote the differential operators soSso tdt respectively. To prove (i), we 
need to show that S'^T^D“^Z(2o, sq, 2 ) is bounded by CsqD^^ for all (j,/c,a). This is clear 
when j = k = \a\ = 0 , directly from a pointwise estimate on FI. We prove by induction on 
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j + k+\a\. We find that w := S^T^D’^^Z{zo, so,t) has component satisfying an ODE of 
the form 


dwi 
dt 


= E 

j 


m 

dz. 


-Wi 


B, 


where i? is a sum of products of factors, each of which is a b-derivative of the form 
applied to H or Z, and where the total number of derivatives applied to any factor of Z is 
strictly less than j + fc + |q!|. By using an integrating factor, the bound on any b- 

derivative of H, and the inductive assumption for lower-order b-derivatives of Z, we deduce 
a similar bound on S^T’^D'^^Z, completing the proof. 

To prove (ii), we write H = i/sm + He, where Hsm is t times a smooth function, and He 
is conormal of order {^ 1 ,^ 2 )- We write Zsm for the solution to the ODE 


(A.4) = Hsm{Zsm, Zo, So,t). 

dt 

Then Z^-m G tC°° using standard ODE theory. So consider Z — Zsen- This satisfies the ODE 

(A.5) ^d{Z^Z^ ^ HsmiZ, Zo,So,t) - HsmiZ^m, Zo,SQ,t) + He{Z, Zo,So,t). 

dt 

It suffices to show that Z — Zsm is conormal of order (/3i,/32). We prove, by induction 
on j + k + |a|, that S^T^D'^^{Z — Zsm) is bounded by When j -|- fc -|- |q!| = 0, 

notice that the RHS of (IA.5I) is bounded by CIZ — Zsm| + We conclude, using an 

integrating factor, that |Z — Zsm| is bounded by CsQ^-t^^. Now consider the b-differential 
operator applied to Z — Zsm. The argument is similar to part (i). Consider the 

ODE satisfied by (Z — Zgm). On the RHS there will be a sum of products of factors 

of various sorts. These terms must of one of the following type. The first type is 

E f ^Hsmi^^Z') Qjrpkr^a ry (-^sm) ry \ 

( sz, ® ^ - 35 —® ^ 

= E - Z.„,) + Z.„,. 


Notice that the first term is a bounded multiple of D^^{Zj — Zsmj), while the second 
is bounded in magnitude by ClZ — Zsm|, and hence by Cs^^t^"^. 

The next type are terms that involve lower-order b-derivatives of Z and Zgm. All such 
terms include a factor that is either of the form (Zj — Zsmj) or {S^ D“^i/sm)(^)— 

(S'-! Hsm){Zsm), or else involve He- Using the inductive assumption, this gives an 

ODE of the form 


dwi 

dt 


E 

j 


dH, 

dzi 


Wj + B, 


for S^T^D‘^^{Z — Zsm), where B is bounded by Csg^t^^. As in part (i), we conclude that 
S^T^D’^^{Z — Zsm) is bounded by 

The proof of part (iii) is similar to part (ii). We write H = Hsm + + Hp, where 

Hsm and H^ are smooth. We first find a function Z(zo,soU) that solves 


(A.6) t— = Hsm{Z,zo,so,t) + sl^r^H^{Z,zo,SQ,t) 

up to an error which is conormal of order (/3i,/32). To do this, we start from the solution 
Zsm of (IA.4I) . and modify it in order to solve away the term Z^{z, Zq, so,t) to leading 
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order, both at sq = 0 and at t = 0. We propose an ansatz of the form Z = Z^m + 
Zj{z, zq, so,t), where Z-y is C°°. Let v{zo,t) be the restriction of Zy to sq = 0i and 
w{zo, Sq) be the restriction to t = 0. To simplify notation, we shall suppress the dependence 
of all quantities on zq from now on. 

To see what the functions v and w must be, we substitute Z^m + Zj(sQ,t) into the 
ODE. This gives a polyhomogeneous expansion both as sq 0 and t —>■ 0, with the first 
possible non-integral power Sq^ as sq —>■ 0 and as t —>■ 0. We seek to make these powers 
agree on the LHS and RHS of the ODE; this will determine v and w uniquely. 

Computing the Sq^ terms of the RHS and LHS of (A. 7) and setting them equal gives 


(A.7) 


^72 I 




a77sm(^sm(0,t),0,t) 
dz, ' 


+ L^7j(-^sm(0, t), 0, t 


Dividing by gives an ODE for Vi which has a smooth solution. Moreover, since iJsm = 
0{t), the the value of at t = 0 is given by 

(A.8 ) t,( 0) =72-'^f7,(^o,0,0). 

Similarly, the coefficient of of the expansion at t = 0 of the ODE is given by 

(A.9) s2;^72W*(so, 0) = 9Hsm(^sn^(5o,0),so,0) ^^ -h77..y^(Zsm(so,0),so,0)y 


Clearly this has a smooth solution Wi{so), with Wi(0) = -(zp, 0, 0) = 'yi(O). Since 

v(0) = w(0), we can find a smooth Zj(so,t) that agrees with z; at sp = 0 and with w at 
t = 0. Then it is easy to check that Zsm + Sg^t'^^Zj(z, zp, sp, t) solves the ODE (IA.6I) up to 
an error term that is conormal of order (/3i,/32), provided that e is sufficiently small. 

To complete the proof, we look for a solution ^'(zp, sp, t) of the ODE 

(Z, Zp,Sp,t) 

of the form Zsm + Zy{z, zq, so,t) + Zp. It suffices to show that Zp is conormal of order 
(/?i,/32 ). This is proved using exactly the same argument as in (ii) above, so we omit the 
details. □ 


In the course of this proof, we have essentially proved the following perturbation result: 

Lemma A.3. Suppose that z,s solve the ODE (IA.2I) . where F,G € tC°°. Let F,G he 
functions in s^^G'^C°° + where and Pi are as in Provosition \A.S\ part (Hi), and 

let F,f = F + F and G» = G -I- G. 

Let Ztf, solve the ODE with F,G replaeed with Ff,G^, and with the same initial con¬ 
ditions as in (IA.2I) . Then 

z(t)-z,(t),—(s-s,) e 
Sp 

The last result we shall need is closely related related to Proposition IA.21 but where the 
initial conditions are specified at a positive value of t, say t = S, where we suppose that 
i5 > 0 is sufficiently small that the solution exists on the time interval t G [0,(5], and we 
are interested in the value at t = 0. To state these results we need to introduce spaces of 
functions with different sorts of regularity in the s and the t variable. We write A^Gf^ 
for the space of functions with conormal regularity of order /3 in the s variable and G°° 
regularity in t and z. 
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Proposition A.4. Let {z, s) solve the ODE 


(A.IO) 


dz 

tL = sG(3,S,() 


^(1) = Zo, 


s(l) = So > 0 


with initial conditions now at t = 1. Then 

(i) Suppose that F,G are as in (i) of Provo sition VKT^ Then 


z(zo,so,t) - zo, —(s(2:o,so,i) - So) e x R+ x M+). 

^0 

(ii) Suppose that F,G are as in (ii) of Provosition [X7B . Then 

z - zo, f (s - so) G C- + x R+ x K+). 

(Hi) Suppose that F,G are as in (Hi) of Proposition\X^ Then 

z - Zo, -(s - so) G + s 2 ^ {G^ + G^G°-) + All fe + + A^C^) • 

The proof is essentially identical to that of Proposition lA.2[ and so is omitted. The only 
difference is that, instead of integrating from t = 0, we integrate from t = 1, so that, for 
example, when we integrate a term of the form with respect to dt/t, we only get 

Cf° + rather than just G^Gf°, accounting for the extra terms in Proposition IA.4I 

compared to Proposition lA.2l 

A simple consequence of Proposition IA.4I is 


Corollary A.5. In case (Hi) of Proposition the functions z(0),s(0)/so are (7°° + 

Sg^C°°’'^ functions of the initial data (zo,So). 

A.2. The sojourn relation. As we describe concretely in the following subsection, accord¬ 
ing to [ 12 ], the Poisson operator is a microlocal object associated to the ‘sojourn relation’ 
SR. We now proceed to describe SR and determine its regularity properties. In this paper, 
we shall take the viewpoint that SR is a Lagrangian submanifold of T*R‘^ x T*S‘^~^, that 
extends nicely to a certain compactification of this spacc0. 

First we describe this (partial) compactificatior0 of r*R‘’* x T*S‘^~^. Let r denote the 
radial variable |a;| and let y be local coordinates on We write h^^{y) for the (dual) 

metric on with respect to these local coordinates. If we write {X,ri) for cotangent 

variables dual to (r, y) on W^, then it is natural to use (A,y = rj/r) near infinity, as these 
are variables that are homogeneous of degree zero under dilations, i.e. remain of fixed length 
as r — >■ oo. We write rj' for a cotangent variable dual to y' on T*S‘^~^, and scale it in the 
same way as 77 ; that is, let p,' = y'/r. Finally, we radially compactify Euclidean space by 
introducing p = r~^ and adding a boundary at p = 0. (However, the space is still not 
compact as /i, p' vary in and </> varies in R.) 


(A.II) We denote the space with coordinates (p, y, y'; A, p, p', cf) by X. 


^In 1121 . the sojourn relation was viewed as a Legendre submanifold of a space with one extra dimension, 
with the extra coordinate being the variable denoted 0 below. Here, we take the view that 0 is a function 
defined on SR. 

^Our partial compactification serves to make the energy surface + V = compact, which is all 

that matters. 
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A more invariant description of this space is given in m, but we wish to avoid the geometric 
intricacies here. The space A is a manifold with boundary, with the boundary defined by 
{p = 0}, and we ignore the apparent singularity in p at r = 0 as we work in a neighborhood 
of p = 0. 

The space X (or at least its interior) is a symplectic manifold with contact form A 
dxj + dp' A dyi- Let V be the Hamilton vector field for the Hamiltonian + V{r,y) — E = 
y? + IpP + V — E, and let V' = rV. In the coordinates (r, p, A, p, t/', p'), V is given by 




(A.12) 


dr 


,’2 dyi 
/9/ibpjpj 91/\ d 

V dyk r2 dyk) dpk ’ 


dr / 9A 


where we sum over repeated indices. In the coordinates (p, y, A, p, y', p'), V' is given by 


V' = -2X( 


(A.I3) 


^ ^ ^ D— 

V dp dp dp' / dfji 

\ 9 fdh'y 

+ (2h p,p^+p )_-( 


dyj 

_ dV\ d 


dyk dyk J dpk ' 


We now perform the operation of ‘blowing up’ X at the submanifold Z = {p = 0, p = 
0,p' = 0}. This operation consists of replacing Z with its inward pointing spherical normal 
bundle, which turns the space X into a manifold with codimension 2 corners that we shall 
denote [X]Z]. This means essentially that \X\ Z] is the ‘minimal’ manifold with corners on 
which the polar coordinates 


Pz — (p^ + \p\'^ + 


PB = —, d = ( —,—), 


Pz 


Pz Pz 


together with p, y', extend smoothly up to all boundary faces. It can be viewed as the 
geometric realization of polar coordinates at Z, that is, the space on which polar coordinates 
are smooth. The space [X; Z] has two boundary hypersurfaces. One boundary hypersurface 
is the original boundary p = 0, or rather the lift of this to the blown up space; we shall 
denote this B. The other is the boundary hypersurface Z created by blowup. We let p^ 
denote any boundary defining function for Z] the above formula for p^ is just an example, 
as any p^ satisfying the properties for bdf’s (see Section rA.l|) will work, and in fact when 
convenient we will take p^ = \p\ near the intersection of B with Z (the ‘corner’) and p^ = p 
in the interior of Z. It follows that ps '■= p/Pz is ^ boundary defining function of B. Away 
from B, coordinates near Z are 


(A.I4) 


p, y, y\ A, 


or equivalently one can take (p^, p, p', p, p', A). Indeed, both p/p = r/ and p'/p = p' 
are bounded maps on compact subsets of the interior of Z, and thus, not only can we take 
p^ = p but the above functions can be checked to yield a coordinate patch on a tubular 
neighborhood Z° x [0,e)p^. 

We now write V' on the space [X\Z]. Notice that pdp + pdp + p'd^/ is precisely Pz^p^- 
Also, we can easily check that 9^^., d\ and p^dfi smooth vector fields on [X] Z], Using 

the assumption that V € p°‘C°°’-{R‘^) C {p 2 ;Pb)°‘C°°'-{\X] Z])^ where we use the notation 
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of Section Ell we compute that V' takes the form 

V' = ( - 2Ap^ + + 0(p“+Vz“-'C“'^)ap3 

(A-15) , gy. 

+ p-^)dx + o{p,c°° + 

where the conormal coefficients of dp^ and dp^ come from the dy^V coefficients of dp^ in 
(IA.13I) . Note that V' is a ‘conormal b-vector held,’ meaning it has conormal regularity and 
is tangent to both boundary hypersurfaces B and Z\ this can be seen directly by noting 
that all the , resp. dp^ terms vanish at B, resp. Z. 

We will multiply V' by a function so that near Z we can use as a parameter for the 
how. Thus for 0 < c < 1/2 to be chosen below, denote by k the function on SR equal 
to 1 for |A| < cx/E, and equal to —2{sgn\)\fEp^/{y'p^) for |A| > (1 — c)x/E. Letting 
V" = Kp 2 ~^V', we have that V''{p 2 ) = —2(sgnA)\^ near U±9± SR, and it follows that for 
small enough c, V” is a smooth vector held on the interior of the blown up space taking the 
form 


(A.16) 


V" = -2(sgnA)V:Eap^ 


in the region |A| > (1 — c)Ve, while the coefficient of dp^ lies in C°° + outside 

this region. Notice that k = 1 at the boundary of SR, and that V" is tangent to B, but 
transverse to Z, pointing ‘inward’ for A < 0 and ‘outward’ for A > 0. We also note for future 
reference, that, if (V°)" is the corresponding vector held for the zero potential, that 

^ V" - (V°)" = 0{p%p^^-^C^^^-)dp, + 0{p%p^^-^C^^^-)pBdp, 


Definition A. 6. We define the Lagrangian submanifold SR as follows: we start from the 
‘initial eondition’ 

(A.18) d- SR := {p^ = 0,y = y\p = -r?',A = -Ve} c [X;Z], 

written using the coordinates (IA.14I) . which is a submanifold of zE Then SR is defined as 
the flow out from d- SR using the vector field V", that is, the union of all integral curves of 
V" starting at points of d- SR. 

Lemma A.7. All integral curves ofV" starting at d- SR reach the set Z n{X = +^/E} in 
finite time. 

Definition A. 8. We define d+ SR to be the intersection of SK with Z (1 {X = y/E}- 

Remark A.9. In |12j the sojourn relation was described as having conic singularities at the 
outgoing radial set G'^, which were resolved by blowing up the span of this set. This blowup 
corresponds to the blowup of Z already performed here. 

Proof of Lemma [X^ Notice that d- SR is contained in the energy surface {A^ + pipj + 
V = E}. By conservation of energy, the integral curves starting from d- SR are completely 


®Near the boundary of Z, we use the coordinates llA.llll and write it in the form 
(A.19) {|/r| = 0, —p' = p, |/r//|/r| = 1, oj = uj', A = —VE, <j> = 0}. 

This is clearly a smooth submanifold of Z. 
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contained in this energy surface. Noting that ^ = 0 and V = 0 at Z, the integral curves can 
only meet Z at A = iL\/E. 

We first show that trajectories contained in the original boundary hypersurface B return 
to Z (1 {X = +'/E} in finite time. In this region, since p/\fJ,\ and |Af^|/|/i| are bounded, we 
can take the boundary defining function for Z to be = \p\- Then consider the vector 
field P 2 ~^V', which is the same as V" up to reparametrization, and hence has the same 
integral curves. We compute that inside the boundary hypersurface B, the variables A and 
1^1 satisfy 

(A.20) a =1^1, 1^1=-A, |^'| = -A, ImI = 

This has an exact solution \p\ = \p'\ = Vasins, A = —y/Ecoss where s £ [0,7r] is the 
‘time’ parameter along this reparametrized bicharacteristic and exp is the exponential map 
on the sphereQ In particular, it returns to Z in finite time, at A = +y/E, as claimed. Then 
by continuity, nearby trajectories also reach Z in finite time. As observed above, this can 
only be at A = Ey/E and by continuity, it must be at A = +y/E. 

Now consider integral curves starting at d- SR and in the interior of Z. These integral 
curves immediate pass into the interior of [T’;Z], i.e. into {p > 0}. By the nontrapping 
hypothesis, they return to {p = 0}, and this can only be at Z, as the vector field V" is 
tangent to B. Since V" is inward pointing at Z for A < 0 and outward pointing for A > 0, 
according to (IA.16I1 . this must occur at A > 0, hence at A = +y/E. □ 

Note that the interior ot Z D {X = y/E} can be identified with x indeed 

it as discussed above, (y, r], y\ rj') give smooth functions on the interior of p^ = 0, and thus 
ZC\{X = yfE} inherits a symplectic structure. The interior of the set d- SR is thus identified 
with as the diagonal, and d- SR itself is in fact the ball bundle obtained by radially 

compactifying the fibers of The boundary d+ SR of SR restricts to be a Lagrangian 

submanifold of this space, and coincides with the graph of the reduced scattering map Se 
at energy E. Thus, the submanifold 9+ SR is precisely the Lagrangian for the ‘absolute 
scattering matrix’. The relative scattering matrix Sh, which is the object we are studying 
in this paper, is the composition of the absolute scattering matrix with the antipodal map 
multiplied by this normalization ensures that the scattering matrix for the zero 

potential is the identity. 

Note further that the integral curves of V" have initial condition on the manifold with 
boundary d- SR, that the boundary of 9_ SR can be defined by the restriction of ps to 
d- SR (as its boundary is exactly its intersection with i3), and that one expects integral 
curves 7 p(t), where p € d- SR is the initial value, to not be smooth in three places: 1) at 
T = 0,2)atp£i3, and 3) when r = Tp, the exit time., i.e. the time when 7p intersects 
a+SR. 

We can now prove Proposition 12.II 

Proof of Provosition \2. 1[ Let T{y', rj') be the time (in terms of the vector field V") taken to 
reach d+ SR starting at (j/', p') £ d- SR. Also, let Y (y', y', r) and N{y', p', r) be the solutions 
of the ODE (IA.I6|) for y, respectively p. For large |y| we use inverse polar coordinates 
p, \p\~^ and write N and for the corresponding ODE solutions. Thus, the map S can 

be expressed in the form 

(A.2I) 5(y', p') = {Y{y', p' , T{y', p')),N{y', rf, T{y', p'))) . 


^As this is happening, y' traces out a geodesic on of length tt, i.e. half a great circle. 
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We first prove the following claim: 

(A.22) 

For \ri'\ < R < oo, Y{y',r]',T{y',r]')) and N{y',rj',T{y',r]')) are C°° functions of {y',r]'). 
For large \ri'\, (Y, N) are C°° + functions of [y',f]', 

while |7V|-i is \r]\-^ times a C7°° + function of {y',fi', 

To prove this for \r]'\ < R, we choose a small <5 > 0 and write T'{y',r]') = T{y',T]') — S. 
Because of the form of V” near d+ SR, this is the time taken for the trajectory starting at 
(?/^ v') G SR to reach the set {p^ = S,X> 0}. As a consequence of ProDOsition lA.2l we see 
that Y{y',ri',T'{y',jf)) and N{y',r]',T'{y 'are C°° functions of {y',jf). (Unfortunately, 
we cannot immediately make the same claim with T replacing T', because Y{y',ri'^ t) fails 
to be smooth in r precisely at r = T.) Now define the map 

(?/o,?7o) Ss{yo,'no), 

where Jyg]rio trajectory that meets i9-|-SR at (yo^Vo) and SsiyoTVo) are the (y,??) 

coordinates of the intersection of ^dh {p^ = S,X> 0}. Again using Proposition IA.2[ 
we see that Ss is a smooth map. Moreover, since V" is Lipschitz, the map is invertible 
for 8 sufficiently small. From these observations, and (IA.21I) . we see that 

5(y', r!) = (Y{y', y', T'{y', y')), N{y', p', T'{y', p'))) 

is smooth. 

To prove the claim for \p'\ large, we follow exactly the same steps, replacing C°° regularity 
by C°° + regularity in terms of the boundary defining function \p'\~^ for d- SR, 

making use of Lemma lA.ll 

Now we consider the effect of the potential V (compared to the zero potential) on these 
functions. Let A be the antipodal map on the sphere, and A* the induced map on its 
cotangent bundle. In the case of zero potential, at d+ SR, A*{y, p) is equal to (y', p'). Since 
the potential V has the effect of perturbing V" by an 0{p%) term, pB = \'ri'\~^, we see 
from Lemma [A.31 that y is given by A{y'), where A is the antipodal map, plus a \p'\~°‘C°°'’^ 
function of the initial values {y',p',l/\p'\). Similarly, after applying A*, l/\p\ is equal to 
1/ly'l plus a \p'\~°‘~^C°°’'^ function of {y\p',l/\p'\). If we write these statements in terms 
of the Euclidean variables p and p', they translate precisely into (12.61) . 

The function p is discussed in (IA.30I) - (IA.31I) below. 

□ 


A.3. Semiclassical parametrix for the Poisson operator. Heuristically speaking, e.g. 
from (na), the Scattering matrix is the limit of the incoming Poisson operator to the sphere 
at infinity after suitably rescaling, composing with the antipodal map, and localizing in 
frequency so as to extract only the outgoing part. We make a more precise statement now 
followed by a characterization of the Schwartz kernels of both the Poisson operator and 
scattering matrices. 

The Schwartz kernel of the scattering matrix, as a half-density, is the distributional limit 


of 




where Ph = Ph{E) is the incoming Poisson operator, A is the antipodal map, and Mout is 
a cutoff to semiclassically outgoing frequencies. This requires some explanation, which we 
give a rough version of now with details to follow. First of all, here we are regarding Ph has 
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a half-density, by multiplying by = \r'^~^drdujdu}'\^^‘^ on x For example, 

Ph.oiE), the incoming Poisson operator for the zero potential is 

Ph^oiE) = 

By [Ml Eqn 1.13], for ip € letting A denote the antipodal map of we have 

Ph,o(E)(iPldujl^''^) ~ + 

and for general (decaying, smooth) potentials, Ph{E){ip\dui\^/‘^) satisfies the same expression 
with A*Ip replaced by A*Ship, where Sh = Sh{E) is the scattering matrix. Taking into 
account the half density factor on one then has 

r-i/2e-*’''/®/'‘P,, (E) (V'Idcu 11/2) 

_ (e-2*’'%/®Ae3"('^-iV(a;) -7e-3’^*(‘^-iU*S',,(i/-))| —dw|i/2. 

r 

The half-density |c?r/r|i/^ is special; it is exactly the radial half-density which makes sense 
to leading order invariantly at the sphere at infinity of the radially compactified Euclidean 
space One thus wishes to cancel off the |c?r/r|i/^ factor, to microlocalize away from the 
^- 2 irVE/h frequency, and then take the limit r —>■ oo. Composing with will then 

give the scattering matrix. 

In the case that E is a smooth function viewed on K'l, which requires in particular 
that a is an integer, the semiclassical Poisson operator was constructed in [12] as a sort 
of ‘boundary value’ of the resolvent kernel. However, the Poisson operator can also be 
constructed directly. Here we make some remarks on this construction in the case that V 
has regularity p“C'°°’C 

We wish to construct a Fourier integral operator F which is a parametrix for the Poisson 
operator. That is, it should have the property that (h^A + V — E)Fh{(p) € h°°p°°C°°(R‘^), 
and also that 

FhiP)) ~ r-(‘‘*-i)/2(e-*^’-/^<)>(w) + F^^/'^iPi-uj)) + o(r(‘'-i)/2), r = \x\ ^ oo. 

That is, up to 0{h°°p°°C°°) errors, Fh<p is a distorted plane wave for A-l-E of energy E, and 
has incoming boundary data (p- Then we will have = Sh{(p) up to an 0{h°°C°°) 

error. 

Based on m, our ansatz is that Fh is a Fourier integral operator associated to the 
Lagrangian submanifold SR. The principal symbol ctq should satisfy the transport equation 

£v«o = 0, 

where Oq is a half-density on SR. Moreover, at 9_ SR, we have an initial condition for 
Oq. This arises from the microlocally incoming condition on the plane wave Fh<p, that is, 
the condition that the incoming boundary data be (p. Microlocally this translates to the 
condition that p^^^ao restricts to d- SR to be the canonical half-density \dy'dr]'\^^^ there. It 
is not hard to see that this implies that the half-density ao is equal to \dy'dr]'dt\^^^, where 
t is the time parameter along the V-trajectories. 

Now consider the higher order symbols in any local parametrization of the FIO. A local 
parametrization takes the form 

•/ j=0 P 


(A.23) 
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Without loss of generality, we can assume that aj depends on a minimal number of variables, 
that is, dim SR = 2d — 1 of the variables {y,y',p,v). We denote these variables collectively 
by A. In that case, the symbols aj, j > 0, given by 




d{\,dy-^) 


- 1/2 


MA|V2, 


d{y,y',p,v) 

are formally determined by (Tq and are solutions of an equation of the form 


j — Q{7j — 

where Q is a second order operator on half-densities on SR, depending on the particular 
variables on which aj depends. The operator Q is induced by the Laplacian. Because of 
this it is times a b-differential operator of order 2. The regularity of the coefficients of Q 
on SR is determined by the regularity of SR, i.e. they take the form C°° + p^°‘~^p%C°°. 

We can write this transport equation on SR in the coordinates {y\ri\t). Writing aj = 
Sj\dy'dr]'dt\^^^, it implies the ODE 

d 2 ' 

— Sj -|- ksj = p Qsj-i, 

where Q is a scalar second order b-differential operator and, near d± SR, we have k = 
2p^dpA{y', y', p). Now changing variable from t to dividing by a factor of p, we have 

near d± SR, 

(A.24) 2Kp2dp^Sj = {2pdpK(ii', p))sj + pQsj-i. 

Moreover, we have an initial condition Sj = 0 at 9_ SR. Propositions IA.2I and IA.4I apply 
to this ODE and show that Sj has the regularity + Pz°‘~^P b^°° away from d+ SR, 
and has the regularity given by part (iii) of Proposition IA.4I near d+ SR. In fact we can 
conclude more regarding the vanishing of the Sj using the ODE comparison lemma. Lemma 
IA.3l above. Indeed, as the coordinates {y', y', t) provide coordinates on both SR and SR°, the 
free scattering relation, we can compare the Sj with the s®, the functions arising analogously 
in the free case, which satisfy that = 0 for j > 1 and Sq = I. It is straightforward to 
check that the ODEs for the differ in an Pz°‘~^P b^°°’^ manner, and thus the lemma 
implies that the Sj for j > I and 1 — so lie in Pz°‘~^ p%C°°^ 

Using these symbols we can build an EIO parametrix for the Poisson kernel, that solves 
the equation (A -|- U — E)Fh4> = 0 up to an error in 

ditional step reduces the error to jl. The error term can be solved away by 

applying the outgoing resolvent, and this contributes a correction term to the parametrix of 
the form w'). This follows from [171 Section 12] (showing smooth¬ 

ness in the y,y' variables) and |2T] (showing that the correction is 0{h°°)). 

Therefore this contributes a smooth term that is 0{h°°) to the scattering matrix, and 
this has no effect on the conclusion of Lemma 12.31 So, in the proof in the next subsection, 
it suffices to analyze the parametrix for the Poisson operator. 


®Further analysis shows that in fact Sj vanishes to order a + j for j > 1, but we do not need this fact 
here. 

®This is just as in the parametrix construction of Melrose-Zworski 
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A.4. The scattering matrix: proof of Lemma 12.31 The scattering matrix is obtained 
by taking a distributional limit of the Poisson kernel as r —?► oo. We can break up the 
Poisson operator microlocally into a piece microsupported away from 9+ SR, and a piece 
microsupported away from 9_ SR. For the first piece, localized away from 9+ SR, if we mul¬ 
tiply the kernel by and then take the canonical restriction to p = 0, we obtain 

the identity operator; that is, if we let this piece of Ph{E) operate on a smooth function 
/(w'), then multiply the result by and then take the canonical restriction to 

p = 0, we obtain /. For the second piece, localized away from 9_ SR, if we multiply the 
kernel by pi/2e“*^/(^p) and then take the canonical restriction to p = 0, we obtain the 
scattering matrix. 

We now see how this happens at the level of kernels. We need to analyze a microlocal 
representation of Fh more carefully near the corner of SR at the intersection of 9+ SR and 
Pb =0. In this region, the Poisson operator can be expressed as an oscillatory integral using 
the results of [12] ■ After a rotation of coordinates, we can assume that pi and 1/pi furnish 
local boundary defining functions p^ and ps, respectively. We can then use coordinates 

(A.25) Z = (2/',pi,l/pi,pj/77i), 2<j<d-l, 

on SR near a point at i? fl Z, that is at ps = p^ = 0. In terms of these we can write the 
other coordinates as smooth functions: 


(A.26) 


yj=Y,{Z), $ = cl>(Z). 


This $ is the same as that in the previous paragraph thought of as a function on SR. In the 
expressions below we will replace pi by a and Pj/pi by n = {v 2 , ■ ■ ■, Vd-i)- Then, according 
to [lH Section 6.3], SR has a local parametrization of the form 
(A.27) 

d-l 

'^{'r,y,y',a,v) = $(?/',cr, -,v) + -(pi - Yi(p',cr, -,v)) + V -Vj{yj - Yj{y',a,-,v)), 

<j p a ^' P O’ 


Then a microlocal parametrix for the Poisson operator Ph{E) takes the form 


(A.28) {VE/{2^h))-^<^-^y^ j 

where ik is as in (IA.27lP^ and the amplitude a has regularity as determined by the regularity 
of the functions Sj as described above. 

If we multiply (IA.28P by p^/'^e~'''^Eph) and restrict to p = 0, using dcrdi; = 
we obtain the following oscillatory integral expression, with the given subsitution we obtain 


(A.29) 






Regarding $ as a function on SR, its behaviour as we approach 9+ SR was worked out in 
[SI Section 2] (where the function is called (j)o)- It is shown that 


$ = VEr + y} + o(l), Pz 0, 


^^The powers of p and a are as given by |12l Equation (6.21)], after taking into account that the half¬ 
density used there is times the half-density \dxduj'\^^'^ used here, times 
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where is a function on SR constant on trajectories, given by 

(A.30) Pb) = X ■ VR, for 7 the trajectory with initial condition (y', f)', pb)- 

Thus ip is the boundary value of a function ip on SR satisfying the ODE V{(p) = x ■ W 
with a zero initial condition at 9+ SR. Using Propositions IA.21 IA.4I and Corollary lA.Si we 
see that 


(A.31) pG p‘^-^C°°’‘^{d+SR). 

(This is the final part of the proof of Lemma EdT ) 

Multiplying by and taking the distributional limit at p = 0 then 

gives 

(A.32) (/E/(27r/i))-('^-i)/2y' ^’^a{0,y',p, h) dr]\dydy'\^/^. 


We write the phase in this equation as (y — A{y')) ■ 77 + G{y',r]). From (IA.3ip . as well as 
Y = A{y') + 0(p“C'~’'=) from ([211), we see that G e p“-iC'~’C 

Then, from the fact that the principal symbol of the scattering matrix is \dydp\^^'^ (see 
[Sj Lemma 3.1]), we find that 


a( 0 ,y',y, 0 ) 


det 


diy,v,y - A{y') 


d^G) 


diy',y,v) 


1/2 


det(Id+4^G)|^/^ 


This shows that 

a( 0 ,y', 77 , 0 ) = l + 0(p“G°°’^). 

Moreover, it follows from the construction of the functions aj in (IA.23P that each term in 
the expansion of a in /i as /i 0 is bounded by 0{p°‘G°°'’^) and also the error contributes 
to lower order that a — 1 G S~°‘. Composing with the antipodal map gives (12. 8 |) and thus 
Lemma 12.31 follows. 


Appendix B. Powers of the scattering matrix 

In this section we will prove Lemma 12.51 which gives an expression for powers of the 
scattering matrix near ‘fiber infinity’. Indeed, recall that, notation as in Lemma 12.51 
Sh = Fi+ F 2 where Fi is a semiclassical FIO with compact microsupport and F 2 (or rather 
its Schwartz kernel) is given by the oscillatory kernel expression in (12. 8 |) , and as discussed in 
the proof of Lemma l4.11 the only non-compactly microlocally supported term in is F^. 

Recall the discussion in Section |4] (see near (14.311 1 explaining that the oscillatory integral 
giving F 2 is to be thought of as a half-density on x Thus two oscillatory integrals 

(B.l) A(y,y') = (27r/i)-(‘'-B J y', y, h) dy 

with 7 = 1,2 define Schwartz kernels Ii{y,y')\dydy'\^^‘^ whose composition as operators is 
given by 

Iiol 2 = {J Ii{y,y")l 2 {y" ,y')\dy"\)\dydy'\^/'^. 

Since this just amounts to integration in the y" variable in all the expressions below we will 
drop the half density factors. 
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Proof of Lemma WfR Suppose that Gi{y, 77 ) and G 2 {y, y) are symbols of order —j3 for /3 > 0. 
We will show that the composition of two FIOs, with phase function 

(B.2) y', y) := {y - y') ■ y + Gi{y', y) 

is an FIO with phase function {y — y') ■ y + Gi{y',y) + G 2 {y',y) + E{y',y), where if is a 
symbol of order —2/3. Thus let be oscillatory integrals as in (IB.II) . with as above 
and with amplitudes Oi € S'™*. The composition has a representation of the form 

(B.3) Ji o /2 := y ai{y,y'' ,y,h)a 2 {y" ^y' ,y' ,h) dy dy dy'\ 

where 


<&(//, y", //', y, y') = (y - y") • y + G,iy", y) + {y" - y') ■ y' + G 2 {y", y'). 

We eliminate the variables {y",y') up to an 0{h°°) error, by replacing them with their 
stationary values and applying the stationary phase lemma m- This works as the Hessian 
of $ with respect to (y", y') is 


0 Id 
Id 0 


+ 0{\y' 




and is therefore invertible, with uniformly bounded inverse, for large \y'\. The stationary 
points in y",y', i.e. the points where Dy^ = 0 : occur at 

y" = y' - dy,H{y'\y',y), 


y = yPdynB{y'' .rf ,y). 


(B.4) 
where 

H{y\rf,y) = G^{y\y)PG2{y\Ti). 

The second line in the above equation array shows that on the critial set we can write 
y = y{y' tU") with y — y' G S^~^. Thus, we want to invert the transformation 

^-dri>H{y",y'y 


(B.5) 


+ 


dy"H{y",y') 


when \y\ is large. It is easy to see, by the method of successive approximations for example, 
that the inverse exists for large y. We claim that the inverse map, which we write in the form 
y"{y'7 y)j y'{y'j y) j is the identity plus a symbol of order —/3. To see this, we differentiate 
(jB.5[l with respect to y' and y to obtain 


(B. 6 ) 

This shows that 


Id 

0 



-dl,y„H{y",y') dl,y,H{y'',y') 

df 


H{y\y') 


y V 



,H{y",y')J \ ^ ^ 


/^eii(y',r?) ei2{y',y) 

\e 2 i{y',y) e 22 {y',y) 


where repeated differentiation of (IB. 6 |) shows that 

feii{y',y) ei 2 {y',y) 
\e2i{y\y) e22{y\y) 

This proves that we can write 


S-i-^ 

s-d 


S-2-t 

S-l-l: 


(B.7) 


y' 


+ 


( h{y'7 y) 
\f2{y',y) 


hGS-^-dj^GS-d. 
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We now write the function $, the restriction of $ to the critical set {£>,,'$ = 0}, 

.gg. 5 = {y - y') ■ V + Gi{y',v) + G 2 {y',v) + E'{y',v), 

E\y,y',rj) = {y'-y”){r,-rj') + G2{y”,v')-G2{y\v)+G,{y”,y)-G,{y',rj) 
Then writing G 2 {y”,v') - G 2 {y',v) = G 2 {y",v') - G 2 {y'',v) + G 2 {y',v) - G 2 {y',v) as 

(y" -y')- [ {dyG 2 ){y'' + t{y' - y"), y') dt + {y' - y) ■ [ {dyG 2 ){y', y + t{y' - y)) dt, 

Jo Jo 

and 

Gi(y", y) - G,{y', y) = {y" - y') • [\dyG 2 W' + t{y' - y"), v') dt, 

Jo 

we see from (jB.7l) that E' is a symbol of order —2/3, so by stationary phase applied to (IB.3I1 . 

ho I 2 = {2TTh)~G-''^') J ’^')/'^b{y, y', y) dy, 

where b{y,y',y) = ai{y,y'',y)a 2 {y",y',y') restricted to the the y'',y' critical set of and, 
as is standard, b € 5'™^“''’”^ with principal symbol given by the product of the principal 
symbols of oi and 02 - 

Lemma 12.51 now follows by applying the above results to repeated compositions of E 2 
(for fc > 1) or F 2 (for fc < 1). Indeed, note that the lemma is already proven for k = 1 
by (12.81) and for fc = — 1 by (I2.15p . We focus on the k > 1 case as the fc < 1 case is 
completely analogous. Assuming by induction that Lemma l^TSl folds for E 2 , consider = 
E 2 o F 2 . Thus the Schwartz kernels of these operators are given by oscillatory intergrals 
/i ,/2 corresponding F 2 and F 2 , respectively, as in (IB.II) with $1 and Gi corresponding 
to F 2 and $ 2 ,G 2 corresponding to F 2 . Thus Gi{y',y) = kG{y',y) + Ek{y',y) and G 2 = 
G{y,y) + E'{y,y), where G comes from the original phase function of Sh, i.e. it is as in 
(12.81) and Ek,E' £ gi-a-e gQj^^g e > 0. Here as in the arguments above we have 
used Remark [SH to switch the roles of y' and y in the perturbation term of the phase 
function. Thus the above arguments imply that the composition has phase function $ = 
(fc + l)G(y',77) + Ek{y',y) + E'{y',y) + E''{y',y) where E" is a symbol of order 2(1 — a). 
Also, the amplitudes satisfy ai — 1 £ 5^““ and 02 — 1 £ then 6 — 1 £ S^~°‘ as well. 

Indeed, this follows immediately from writing y", y' in terms of y' ,y using (IbTI) . □ 
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